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Abstract 

We present a comprehensive model to account for the behavior of suspended 
nanoparticles under magnetohydrodynamic conditions. The Lorentz force 
not only drags nanoparticle floes toward the walls reducing the distance be- 
tween floes resulting a more negative total pair interaction potential energy, 
but also produces extra magnetic-induced stresses inside a floe leading to a 
change of pair interaction distance thus giving rise to a less negative total po- 
tential energy. The model explains quite well the recent experimental results 
showing that magnetic field assists aggregation in laminar or weak turbulent 
flows, but favors floe disruption in turbulent regime. 

Keywords: magnetohydrodynamic(MHD), laminar, turbulent, floe, 
aggregation, disruption 



1. Introduction 

The stability of colloidal suspensions has drawn a lot of interest in the lat- 
est century. According to the well-known DLVO (Derjaguin-Landau-Verwey- 
Overbeek) theory [H, 0], this stability is governed by the repulsive interaction 
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due to the overlap of electrical double layers, competing with the Van der 
Waals attractive interaction between the primary particles. If the particles 
are transported by molecular motion or mechanical mixing, further floe (loose 
structure that composed by primary particles and fluid) aggregation occurs 
due to the particle collisions jij. On the other hand, during the transport 
process the floes are subjected to unequal shearing forces, especially under 
turbulent conditions, giving rise to surface erosion j3] or floe disruption (if. 
Once the rate of floe breakup becomes equal to the rate of floe growth the 
system is in equilibrium and a steady-state floe size distribution is reached. 

Recent decades, magnetic field has been considered by some research 
groups as another factor enhancing dispersion. Svoboda et al. [6] and Parker 
et al.[7] have mentioned that, in stationary or gently flowing suspensions, 
magnetic fields magnetize the particles yielding attractive interaction be- 
tween magnetic dipoles. The modified DLVO potential curve produces a 
secondary minimum at relatively large interparticle distances, where the col- 
loid becomes unstable and flocculates as loose aggregates. This magnetic 
attraction is, however, too small to be observed. Busch et al. [8( have argued 
that combining a laminar flow with a magnetic field, i.e., under magneto- 
hydrodynamic (MHD) conditions, a much more pronounced effect can be 
obtained experimentally. Based on a hypothesis that the flowing fluid is 
conducting, the authors ascribed this effect to the formation of a boundary 
layer near the walls leading to a much stronger velocity gradient to increase 
the collision frequency and, thus, to increase the average floe size Jsj. Un- 
fortunately, this promising concept has not received its matching attractive 
case for turbulent flow since, to the best of our knowledge, the existing mod- 
els cannot explain adequately the available experimental data. Sometimes 
these models are even in contradiction with the experiment. Therefore, it is 
necessary to develop a convincing model to explain the recent experimental 
results, especially the disruption mechanism of floes under turbulent MHD 
conditions. 

In this paper we present a comprehensive model to describe the magnetic 
effects on a suspension under hydrodynamic conditions. We assume that a 
magnetic field produces magnetic (Lorentz) forces acting on primary-particles 
and hence floes, which changes the distances between the floes and the dis- 
tances between primary-particles inside a floe, thus giving rise to a growth or 
reduction of the floe size. Our model predicts that under MHD conditions, 
the size increases in laminar or weak turbulent flows and decreases in turbu- 
lent regime, which explain quite well the experimental data available by now 
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2. Models and discussions 

As we know from Refs. @, 0, @, 0, 10, when a transverse external magnetic 
field (B x ) is applied, the total energy of interaction of colloidal particles is 

V T = V EL + V VDW + V M , (1) 

where Vel = 4:7ieaipQexp[—Ka(s — 2)]/s is the simple approximate form of 
the repulsive energy of electrical double layer, whereby e is the absolute 
dielectric constant of the medium, a is the radius of interacting spheres, iJj is 
the potential at the surface of the particles, k is the Debye-Hiichel parameter, 
and s is the ratio of the distance of the centers of two sphere Ro to there 
radius a; Vvdw = — j + 7s + ^-^r] is the Van der Waals attractive 
energy between two spherical particles whereby A is the Hamaker constant; 
and Vm — — 327r 9 ° o ^ 3 Bz is the magnetic attraction between magnetic dipoles, 
whereby x is the volume magnetic susceptibility of particles and /j,q is the 
magnetic permeability. 

However, under the same magnetic environment, if the suspension is hy- 
drodynamically transported in a pipe-line system, floes flowing along a pipe 
(x-direction) are subjected to Lorentz forces. As it can be seen from Fig. 
1(a), for the floes with positive charges, the Lorentz forces always drag the 
floes toward the walls. This process is similar to the sedimentation effect 
that can be called as floc-migration (FM). 



On the other hand, floe rotation in laminar flows [U] and velocity fluc- 
tuation in turbulent flows are inevitable. In most cases, the former 
one is neglected while the latter should be taken into account. As mentioned 
in Ref. |5|, random velocity fluctuations may induce pressure difference on 
opposite sides of a floe, which causes its deformation and eventual rupture 
[5|, i.e., it increases the distance between the primary-particles inside the 
floe. Thus, as shown in Fig. 1(c), we can expect that in a turbulent flow a 
floe moves along the circular orbit and spins around its center, showing floc- 
s elf-rotations (FSR). This is somewhat similar to the Earth orbital motion 
around the Sun combined with the Earth's spinning around its axis. Such 
FSR causes the constituting primary-particles in the floe to have different 
net vector velocities. In this case we divide the floe into two parts, a and 
b, which have opposite velocity directions u a and w&. When a magnetic field 
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Figure 1: (Color online) (a) Migration of floes towards the walls in a laminar flow, 
(b) Schematics of the distance change between two floes after magnetic field is applied. 
Position- 1 and 2 are the positions of the two floes without a magnetic field. Position- 
1' and 2' are the positions of the two floes with a magnetic field after time t. Obviously, 
distance-1'2' is somewhat shorter than distance-12. (c) Self-rotation of a floe in a turbulent 
flow. The dashed line is a vorticity orbit while the solid line separates the primary-particles 
into two groups, a and 6, which have opposite net vector velocity u a and Ub with respect 
to the applied magnetic field, producing Lorentz forces with opposite directions fi, a and 
fLb- (d) Schematics of the distance change between two primary-particles in a floe after 
a magnetic field is applied. Before disruption, the distance is i?o without a magnetic field 
but ai?o with a magnetic field. It is clear that a > 1. 

is applied, the two velocities produce opposite Lorentz forces and induce an 
extra pressure difference (we identify it as magnet-induced stress) on opposite 
sides of the floe, which enhances the total stress (i.e. the total interaction 
energy Vr between two primary-particles) under turbulent conditions, and 
thus stretches the floe to a more prolate spheroid. Once the total stress 
surpasses the yield stress of the floe, it breaks up. 

In principle, we cannot directly estimate the FM and FSR effects by a 
special potential term as Vel, Vvdw, and Vm shown in Refs. 0,0 because 
the Lorentz forces acted on different charged particles are independent from 
each other. But, by means of the distance change between the floes due to 
the FM effect, or between the primary-particles inside a floe due to the FSR 
effect we can estimate the total energy change in Eq. (|Tj). To do this, we 
need to know the velocity distribution in each position at every time. The 
most precise way is to solve the Navier-Stokes equations by computational 
fluid dynamics (CFD) softwares. But this is a quite complicated task and it 
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goes beyond the scope of this work. 

Fortunately, we can have two simple alternative ways to estimate the 
velocity distributions and hence the distance changes due to the magnetic 
effects. First, for the distance change between floes caused by the FM effect, 
considering that this FM effect is dominant in laminar flows (also possible in 
weak turbulent flows), one can approximate the velocity in the x-direction 
along a 2-dimension (2D) pipe as 

v x = v max (l-4y 2 /D 2 ), (2) 

where v max is the maximum linear velocity at the center part of the pipe, 
and D is the diameter of the pipe. Then, from the physical viewpoint, the 
dynamics of an isolated floe in y-direction is governed by the Newtonian 
equation as 

Fcravity,y FLorentz,y Fviscosity,y TTlf ^ , (3) 

where F Gravity , y is the gravity force, F Lorentz ^ is the Lorentz force, F viscosity ^ y 
is the friction force of the fluid, and mj and y are the mass and translation 
of the floe, respectively. 

If Fiorentz^y « Foravity,y, the system is reduced to the sedimentation case 
and the velocity and moving distance in y-direction, v y and y, are mainly 
dependent on the gravity force. Thus, at steady state, one can easily find 

that [7], V y = 4jL = ^(Ppartide- P fluid) Of, and y = yo + ^{p pa Hicle- P fluid) a}t, 

where p par ticie and pf\ U id are the densities of the floe and the fluid, respectively, 
a/ is the radius of the floe, g is the gravity constant, /i is the viscosity of the 
fluid, and t is the integral time. 

However, in modern nano-dispersion systems the relation F^ orentz ^ y » 
FGravity, y is often satisfied since the gravity effect can be neglected. Similarly, 
by using Eqs. ($2§ and (j3J) the v y and y can be estimated as, 

_ dy _ AD[—1 + tanh(At) 2 ](D 2 - 4yg) 
Vy ~~dt~ 2[D - 2y tanh(At)} 2 ' ( &) 

2y D - DHanh(At) 
V ~ 2D- Ay tanh(At) ' ^ ' 



interactions from Electrical double layer, Van der Waals, and magnetic dipoles are 
not taken into account here since they are balance each other. 
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where A = qNB z v max / (3naffiD), q is the elementary charge constant, and N 
is the total charge number in the floe which is associated with the potential 
of the surface charges, ipo- 

Consequently, combining Eq. fl2]) and (TjJ we can estimate the position of 
a floe flowing in a laminar flow at any time due to the FM effect. Fig. 1(b) is 
the schematic presentation of two floes migrating in a 2D pipe. Position- 1 
and 2 are the originals of the two floes. As it can be seen, without magnetic 
field the y coordinates are not changed, while with magnetic field the floes 
move to position- 1' and 2', respectively, after time t. Obviously, distance-1'2' 
is somewhat shorter than distance-12. To trace in a better way the distance 
change by the magnetic field we introduce a ratio a = dt J^ance-i2 • This a is 
clearly dependent on magnetic field, and potential of the surface charge. For 
such FM case a is smaller than 1. 

Secondly, the distance change between the primary-particles inside a floe, 
due to the FSR effect which dominates in turbulent flows, can also be esti- 
mated in a simple way. Notice that before a floe breaks the extra stress (force), 
for instance, the magnetic-induced stress inside the floe is always relaxed by 
the distance change. For simplicity we only discuss pair interaction between 
any two primary-particles as the classic DLVO theory did to estimate the 
distance change between the pairs. Fig. 1(d) is the schematic presentation 
of the distance change between two primary-particles inside a floe induced 
by magnetic field. We assume that these two primary-particles are balance 
at positions 1 and 2 without magnetic field and at positions 1' and 2' with 
magnetic field. The distances between them are Rq and aRo, respectively. 
For a given turbulent dispersion system without magnetic field, if we assume 
that the length and time scales are the order of Kolmogorov length and time 
scales, respectively, the distribution of the intensity of turbulence over the 
range of wavenumbers (or over the range of eddy sizes) can be determined 



as follow 12 
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u\ 2 = E(k,t)dk, (5) 



where Ui ? is the fluctuation velocity, k is the wavenumber (reciprocal of the 
characteristic length r), and E(k, t) = const, x e 2 / 3 k~ 5 / 3 f or n <^ r <ti L while 
rj is the Kolmogorov length scale and L is the large scales. Therefore, before 
disruption the distance between any two primary-particles, Rq, inside a floe 
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is governed by the force balance equation acting at one primary-particle as, 

\m ia u\ 2 dV T 
Ro/2 ~ ds~ 

where my is the mass of primary-particles. The left-hand side is the force 
caused by the kinetic energy of the turbulence while the right-hand side is 
the static force caused by the total interaction potential. Similarly, for a 
system with a magnetic field, the force balance is changed to, 

(7) 

s=aRo /a 

where u' x 2 is the fluctuation velocity with a magnetic field, and Nq is the 
total charge number in the primary- particle. The second term in the left- 
hand side is the Lorentz force. Furthermore, before disruption, the change of 
the total interaction potential must equal to the change of the kinetic energy 
(energy conservation), 

V T (aR ) - V T (R ) = imxK 2 - u\) + ^m 2 (u' 2 2 - u\). (8) 

For simplicity, we further suppose that the primary-particles are identical, 
that is, mi = m-i = m pp , \u±\ = \ui\ = \u\ and \u[\ = \u' 2 \ = \u'\. Thus, 
combining Eqs. (JTJ, (jBj), ((7|), and (jSJ) we can calculate the four unknown 
variables, Vr, Rq, u', and a. Obviously, this a is dependent also on the 
magnetic field, Reynolds number, and potential of the surface charge. In 
this case a is greater than 1, which means the floe is getting more prolate 
and less oscillating before it is broken apart. 

Now we turn to the discussion of the effect of an external transverse 
magnetic field on the pair interaction potential curve. Fig. 2 presents the 
total energy with respect to the distance between two primary-particles. The 
solid line without symbol is the classic DLVO curve without magnetic field. 
Solid curves with open star, triangle, square, and circle are slight distance 
change cases, such as, a = 0.985, 0.995, 1.005, and 1.015, respectively. From 
the figure we can find that in laminar flow (a < 1), due to the FM effect, the 
curve produces a more pronounced secondary minimum which means that the 
interaction energy is more attractive than the one without magnetic field. In 
this case the floes flocculate more easily from a relatively large interparticle 
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s=Rn!a 
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dV, 



T 



ds 
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Secondary minmum points 



2.2 

Distance s=r/a 



Figure 2: (Color online) Total energy with respect to the distance between two primary- 
particles. The solid black line without symbols is the classic DLVO curve without mag- 
netic field. Solid curves with open star, triangle, square, and circle are the cases of 
a = 0.985, 0.995, 1.005, and 1.015, respectively. Black filled dots in each curve are the 
secondary minimum points. 



distance, giving rise to a relatively large floe size, which is in a good agreement 
with the experimental results reported in Refs. 0, 0, 10. While in a turbulent 



flow (a > 1), due to the FSR effect, the curve is completely above the 
classic one which means that the total energy is less attractive than the one 
without magnetic field. In this case the original aggregates are elongated into 
prolate spheroids by the extra magnetic-induced stress and eventually break 
up into smaller floes leading to a relatively small average floe size. This is 
the reason why the size of floes is reduced in strong turbulent flows reported 
m Refs. @, Hoi Furthermore, in weak turbulent flows, both the FM and FSR 
effects are comparable. Thus, the floe size is dependent on the competition 
between FM and FSR. Using this model we can explain why the floe size 
increases in weak turbulent region in Refs. I, Hol : the FM effect in this case 
is stronger than the FSR. 



3. Conclusions 

To conclude, a transverse magnetic field results in Lorentz forces acting on 
primary-particles and hence floes, which drag floes toward the walls and/or 
elongate floes into prolate spheroids. The former mechanism mainly works in 
laminar flows leading to a decrease of the distances between the floes, while 
the latter dominates in turbulent flows thus giving rise to an increase of the 
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distances between the primary-particles inside a floe. This model predicts 
that under MHD conditions, the size increases in laminar or gentle turbulent 
flows and decreases in turbulent regime. We believe that the presented model 
is not only applicable for physics and chemical engineering science, but also 
for the biology, medicine, etc. A good example here is the investigation of 
the stability of protein molecules. It should be noted here that the model 
in this paper is only one eddy for the regime of turbulence. Further stage 
calculations one should take into account many eddies of different sizes and 
a statistical description should be used. 
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